The Possible J PC I G = 2++2+ State X(1600) 
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I. INTRODUCTION 
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The interesting state X(1600) with J PC I G = 2 ++ 2 + can't be a conventional qq meson in the 
quark model. Using a mixed interpolating current with different color configurations, we investigate 
the possible existence of X(1600) in the framework of QCD finite energy sum rules. Our results 
indicate that both the "hidden color" and coupled channel effects may be quite important in the 
multiquark system. We propose several reactions to look for this state. 
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There has been important progress in the experimental search of multiquark hadrons since last year. LEPS col- 
laboration reported evidence of the + pentaquark with S = +1,B = +1 and the minimum quark content uudds 
0. Such a state is clearly beyond the conventional quark model if it is established by future experiments. A recent 
review of pentaquarks can be found in Ref. Q • 

BABAR, CLEO and BELLE collaborations observed two narrow charm-strange mesons D s j(2317), Z? s j(2457) below 
threshold Is. These states are 160 MeV below quark model predictions. Some authors speculate they are four-quark 
states 0,U These states may admit a small portion of DK or D*K continuum contribution in their wave functions. 
I" ■ But the dominant component of £> s j(2317), D s j(2457) should be cs Q- This topic is reviewed in Ref. Q. 

BELLE collaboration discovered a new narrow charmonium-like state X(3872) in the J/i(m + TT~ channel ||- Its 
mass is very close to DD* threshold. Its production rate is comparable to those of other excited charmonium states. 
Q_i| Very recently BELLE collaboration observed the same signal in the J/ift u channel 0- Its mass and decay pattern 
JL . seems to favor its interpretation as a deuteron-like DD* molecule |l0j . Naively, one would expect the lower production 
rate for DD* molecules. 

C"| , This year SELEX collaboration reported a narrow state D s j(2632) above threshold Its dominant decay mode 
is D s ?]. Such an anomalous decay pattern strongly indicates -D s j(2632) is a four quark state in the SU(3)f 15 
. [ representation with the quark content -^^(dsd + sdd+ suu + usu — 2sss)c .12]. Other possible interpretations were 

^ ' discussed in Refs. fl3| . 

In the light meson sector, /o(980)/ao(980) lies 10 MeV below the K + K~ threshold. It's difficult to find a suitable 
position for them within the framework of quark model. So they were postulated to be candidates of kaon molecule 
or four quark states. Recently there has accumulated some evidence of the four-quark interpretation of the low lying 
scalar mesons from lattice QCD calculation 113 • 

The study of multiquark states started nearly three decades ago in the MIT bag model EUS- 

The reaction 

77 ~~ * PP was suggested in the search of qqqq resonances. Later ARGUS collaboration found evidence of a four-quark 

state in the dominant partial wave J P J Z = 2 + 2 in the reaction 77 — > p°p° — > 4tt 01- Now the observed signal was 

named as V(1600 ± 100) with the quantum numbers J PC I G = 2++2+ [lj|. 

In this work, we employ QCD finite energy sum rules (FESR) to explore whether there exists a resonance in the 
jPCjG = 2++2 + channeL 
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II. FORMALISM 

QCD sum rule approach has proven useful in extracting the masses of the ground state hadrons _2Gj . QSR approach 
can yield the absolute mass scale as the lattice QCD formalism. First one starts from the correlation function composed 
of the interpolating current which strongly couples to the hadron which one wants to study. The correlation function 
can be calculated using the operator product expansion (OPE) technique. As one approaches the resonance region 
from large Q 2 , the nonperturbative power corrections become important gradually. One gets the spectral density of 
the correlation function in terms of quark and gluon condensates at the quark level. The hadron mass enters the 
spectral density of the correlation function at the hadron level. With the quark hadron duality assumption, one can 
extract the hadron mass. 

The construction of a suitable interpolating current is crucial. There are only two independent color structures 
for a tetraquark. Let's first focus on a pair of qq. Since 3 C x 3 C = l c + 8 C , there are only two ways to form a color 
singlet tetraquark from two pairs of qq. Both pairs are either in the color-singlet or color-octet state simultaneously. 
We denote the corresponding interpolating currents as rj 1 ,!] 8 . If we focus on the two quarks, 3 C x 3 C = 3 C + 6 C . 
The anti-quark pair must be in either 3 C or 6 C color state in order to get a color-singlet tetraquark. We denote the 
corresponding interpolating currents as r] 3 ,rf . It's important to note that rf , rf are linear combinations of r) 1 ,rf . We 
refer the reader to the appendix for details. 

We use a general interpolating current for X(1600) which is the linear combination of rf and rf . 

ri^(x) = vl v + Yvlu (1) 
= (Y + l)d l (x)^u l (x)d m {x) 1 ,u m (x) + (Y — l)d l {x)^u m (x)d m {xh„u l (x) + (.g^terms) 

where Y — a + bi is a complex number, a and b are real numbers. The decay constant fx for X(1600) is defined as 

(0|ry^(0)|A(1600)) = /xe^ (2) 

where is the polarization tensor of X(1600) meson. 
We consider the following correlation function 

i J d A xe- 1 ^ < Om^^WHO >= A^ipMp 2 ) + ■■■ (3) 

where 

A^jc/sCp) = \ ( ^ a (p)A v0 (p) + A flfi {p)A lJa {p) - ^A tlu (p)A al3 (p)j , (4) 



A^(p) = g^u - V^Vulv 2 ■ (5) 

In Eq. ©, we have kept the unique tensor structure for J PC = 2 ++ mesons: A^ v ^ a p{p). We note in passing that the 
qq meson states with J PC = 2 ++ have been studied in |22,|2]|, where L=l orbital excitation has to be introduced. 
For any of its four Lorentz indices, A^ v ^ a p(p) satisfies: 

<f A^^Cp) = (6) 



K,^P) - . (7) 

The non-resonant p + p + intermediate states do not contribute to this tensor structure. 
The scalar complex function II(p 2 ) satisfies the following dispersion relation 

rV) = / p( ; s) . ds (8) 

J s - p- - te 

where p(s) is the spectral density. For a narrow resonance 

p(s) = f x S(s - M x ) + higher states . (9) 
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FIG. 1: The Feynman diagrams in the calculation of the two-point correlation function in Eq. 10. 

At the quark gluon level the correlation function © reads 

i Jd 4 xe~^{\Y + l\ 2 [Tr[ 7ll Sl n (x) la S mi (-x)] x Tr^S™^^ (-x)] 

-Tr{j^ n (x) la S m 3(-x)j u S m (x) Jp S m (-x)} 
+2(\Y\ 2 - l)[Ti[ 7 ^ m (x) la S mi (-x)} x Tr[ lv S in (x) 7 0S n j(-x)} 
-Tr[ 7 ^ m (x) la S m ^-x) 7l/ S in (x) 7f3 S ni {-x)\ 
+ \Y- l\ 2 [Tr[^S(x)j a S(-x)} x Tr[j v S(x)j S(-x)\ 

-Tr[ lfi S(x)~, a S(-x)~f„S(x)~f,3S(-x)] + {a «-► /?)] (10) 

where iS(x) = (0\T{q(x)q(0)}\0) is the full quark propagator in the coordinate space. The Tr denotes the summation 
of both the color and Lorentz indices. Throughout our calculation, we assume the up and down quarks are massless. 
The first few terms of quark propagator is 

if) a h i \ n 1 fi ab Xab 2 

lSab{x) = + pf^ 4 + ™ n ~ n m + ^m {9 ^ aGq) + ■ ' ' 

The relevant terms that contribute to this correlator are represented pictorially in Figure 1. 

After making Fourier transformation to Tl(x) we arrive at Tl(p 2 ). From the imaginary part of Tl(p 2 ) we extract the 
spectral density p(s). 

P{S) = W^T^ 1 + ^ |F|2)s4 ~ 2" . is . ^ ( 2 3| y | 2 - 21 ( y + Y *) + ^> 2 (9 2 S GG) 

+ 18^ (1 _ gl y lX«9> 2 - 14^2 ( 13 l y ! 2 - 3 ( y + y *) - ^)(qq)(g s qaG q ) (12) 
where we have used the factorization approximation for the high-dimension quark condensates. 

III. FESR AND NUMERICAL ANALYSIS 

In the Borel sum rule (BSR) analysis there are two parameters: the continuum threshold sq and Borel mass Mb- 
FESR contains a single parameter so |24|. The dimension of the interpolating currents of the conventional hadrons 
like rho and nucleon is not high. Both BSR and FESR yield roughly the same results. If the interpolating current is of 
high dimension, the working window of Mb does not exist sometimes. In this case, FESR may have some advantage 
over BSR H|. 

With the spectral density, the nth moment of FESR is defined as 

W(n,s ) = / dss n p(s) (13) 
Jo 

where n > 0. With the quark hadron duality assumption we get the finite energy sum rule 



W(n,s )\ Ha dron = W(n, s )\qcd 



(14) 
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The mass can be obtained as 



M2 = i^^+Mo) 

W(n,s ) 



In principle, one can extract the threshold self-consistently from the requirement that the hadron mass has the 

iM 2 
ds 



least dependence on so, i.e., = 0. However, the threshold value sq extracted this way is not necessarily physical 



sometimes. The weight function of FESR enhances the continuum part even more than the weight function e~ s l MB 
in the Borel sum rule. One must make sure that only the lowest pole contributes to the FESR below sq- Otherwise 
the result will be very misleading. To be more specific, a naive stability region in so is no guarantee of a physically 
reasonable value for so. For example, the FESR with an extracted threshold so « 10 GeV 2 is certainly irrelevant for 
the possible X(1600) state. 

Besides the weak dependence on so, we also require (1) the zeroth moment W(0, so) > and (2) the convergence 
of the operator product expansion (OPE). Higher dimension condensates should be suppressed in order to ensure 
a reliable FESR based on the converging OPE series. This convergence requirement is not easily satisfied if the 
interpolating current is of high dimension. As shown in the appendix, neither rf nor rf leads to a converging FESR 
for a reasonable value of sq. Only with a mixed current which is the linear combination of rf and rf , can we make 
the OPE series of FESR converging. 

We use the following values of condensates: (qq) = -(0.24GeV) 3 , {g 2 s GG) = (0.48 ± 0.14)GeV 4 , (g s qcrGq) = -m 2 , x 
(qq), m§ = (0.8 ± 0.2)GcV 2 . To simply numerical analysis, we introduce variables ei, £2 which are defined as: 

« 2 + ^(l-ei) 
« = ~(l + e a ). 



They satisfy two constraints: ex < 1, |1 + e?\ < 6-^/10(1 — ei). Now the spectral density reads 

P(°) = 2 i2 . 2! . ^ ( 8 ~ 5£ i)* 4 - 2 i5 . 15 . ^ ( 216 ~ 115£ i + 7f2 + 7) s2 (9 2 sGG) 
5 _ 2 1 

+ 18-7r 2 £lS ^g) + 288 ■ 7T 2 ( ' 65ei ~ e2 ^11^9sq^Gq) . (16) 

From Eq. (|16fl it's clear that we can adjust the variables ei, €2 to suppress D = 6, 8 condensates. In fact there exists 
a parameter space of (ei,e 2 ), with which there exists a stable FESR plateau in the Mx vs sq curve. The extracted 
so at the stable plateau is physically reasonable. With this so the zeroth FESR moment is positive, and the OPE is 
convergent. One typical set of these variables is e\ — —0.02, £2 = —2. After we divide each piece in the FESR by the 
perturbative term, the zeroth moment reads 

W(o,* )~i--5 -3 ZT- (17) 

Sq Sq Sq 

The positivity requirement leads to so > 2.3 GeV 2 . 

For s > 2.3 GeV 2 , the OPE also converges well. The variation of M x with s is presented in Figure 2. The 
extracted mass is 

Mx = (1.65 ± 0.15)GeV . (18) 

As can be seen from Figure 2, the value of Mx does not change much if we truncate the OPE at D — 6. This is 
another sign of convergence of our FESR. 



IV. DISCUSSION 

The state X(1600) with J PC I G = 2++2+ has inspired many theoretical papers pfjj . The potential scattering of 
via the a meson exchange was proposed to explain the signal in Ref. [23. In the diquark cluster model [2^|, the 
mass and decay width of V(1600) is studied together with other tetraquarks assuming the 3 C color wave function for 
two quarks. Using the potential model A(1600) mass was estimated to be 1544 MeV recently [2^. 

In this work we have investigated this state using QCD finite energy sum rule. With any single current rf , rf , rf , 
rj e , the operator product expansion of the corresponding correlator does not converge. After we consider the linear 
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FIG. 2: The variation of the mass M with the threshold so (in unit of GeV 2 ) for the current Q. The long-dashed and solid 
curves correspond to FESR when truncated at dimension D = 6 and D — 8 respectively. 

combination of rf and rf , we have derived a converging FESR and observe a resonance signal at 1.65 GeV, which 
is close to the experimentally observed X(1600) state. Our analysis indicates both the "hidden color" and coupled 
channel effects may be important in the multiquark system. 

From charge, angular momentum, parity, isospin, C parity, and G parity conservation, we may obtain the allowed de- 
cay modes of A(1600). G parity requires even number of pions in the final state for pure pion final states. For example, 
A°(1600) — > p°ir°, um + iT~ is forbidden by G parity. C parity forbids A°(1600) — * uj(3tt°). Isospin conservation forbids 
the following decay modes: A°(1600) — > ojtt°,ujuj. The possible modes are A°(1600) — » tt q tt°, 7r + 7r~, Att, p + p~, p°p°. 
Angular momentum conservation requires D-wave decay for the two pion mode. Therefore, two pion decay modes 
may be suppressed compared with S-wave pp mode. 

Filippi et al. reported the possible existence of the J PC I G = ++ 2 + state decaying into tt + tt + in the reaction 
np — ► 7r + 7T + 7r _ |3C|. One the other hand, some constraint has been obtained on the possible resonance in the scalar 
isotensor channel from the phase shift analysis of tttt — > pp — > hit scattering [23|. Constraint on the D-wave 1=2 
resonance from the similar analysis of the same reactions will be very desirable. 

The process 77 — + pp — » A(1600) — > An was suggested and used in the search of A(1600). We suggest the following 
reactions to look for this charming state. 

• J/ip decays at BES, CLEO and T decays at BELLE, BABAR 

Symmetry considerations requires A (1600) is produced together with odd number of pions in J/ip decays. One 
pion mode is forbidden by isospin conservation. So the favorable decay chain is J/ip A(1600) + 3tt — > 
2p + 3tt — > 7ir. There is also some chance in J/ip — > A (1600) + 3tt — ► 2ir + 3ir —> 5ir. 

• Hadron reactions especially charge exchange processes 

Charge exchange processes are useful for the production of A ++ (1600) if such a state exists. Some of them are 
7T+ N -> 71-- + X++(1600) + N^5ir + N,p + N->n + N' + ^++(1600) where N is either a nucleon or nuclei. 
The later process can be studied at CSR facility at Lan Zhou. 

• Anti-proton annihilation on the proton or deuteron targets 

Let's take A(1600) production from anti-proton annihilation on the proton target as an example. Since the 
isospin of pp is either or 1, A(1600) should be accompanied by one or several pions. However, kinematics 
allows only one pion decay mode: p + p — > A(1600) + tt° — > 5tt. So we get C p p = (—^ l pp+ s pp = +,I pp = 1. 

(1) When the pp pair is in the S-wave, L p p — S p p — J pp — 0,P pp — — . Angular momentum and parity 
conservation requires £7^(1600) = 2; (2) When the pp pair is in the P-wave, L pp = S pp = 1, P pp = +, J pp — 1 or 
2. We get £^(1600) = 1> (3) When the pp pair is in the D-wave, L pp = 2, S pp = 0, J pp = 2, P pp = ~. We have 

-^7rX(1600) = 0. 
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V. APPENDIX 

In this appendix we list the currents with different color structure and their spectral densities. 
The color-singlet p + (x)p + (x) type interpolating current reads 

ril v (x) = d\x) ltl u\x)di{x) lv ui(x) - l/ig^d 1 (x) lp u i {x)${x)Y^{x) , (19) 

which "i, j" are the color indices. 

The current with color octet structure reads 

ijjUs) =d\x) (^j 7 ^(x)d m (x) (j) lv u n {x) - \g,J\x) (^j lp v? (*)<?» (j) Yu n {x) (20) 
where 4r is the SU(3) C generator. 



Using the identity of the A matrix 

2 ' 2 2 



we can rewrite Eq. I|20|l as 



rilM = -\d\x) ltl u l (x)d m (x) lv u m {x) + ^d 1 (x)^u m (x)d m (x) 1 ,u l (x) + (pterins) . (22) 

Alternatively, when two quarks are in the 3 C state, we have 

Ti*l(x) = d l (x) 1 ^u l (x)d m (x) lly u m (x) - d}{x) 1 , 1 u m {x)d m {x) lv u\x) + (^terms) . 
The current with the 6 C color structure is 

rfc(x) = ${x) ltl u\x)d m {x) lp u m {x) + ${x) ltl u m {x)d m {x) lv u\x) + (ff^terms) . 
The spectral densities of the above four currents 

As) = - W^^GG) + 2^ S ^ + 28^<«*><*S^> ' (23) 

1 1 97 1 131 

1 23 1 13 

P 6c ( s ) = 2 n .21-7r6 a4 - 2 i4. 15 . 7r6 s2 (g' GG ) " ^2 S ^ 9 > ~ 144 ,^ 2 (qq)(g s qaGq) . (26) 

In order to study the convergence of OPE series in FESR, we divide each piece in the zeroth moments by the four 
quark condensate. 

W x (0, s ) ~ 7.5 x lO^s 3 , - 4.7 x 10~ 2 s + 1 - — . (27) 



W 8 (0, s ) ~ 3.5 x lO^s 3 , + 8.9 x 10~ 3 s + 1 
W S (0, s ) ~ 2.7 x l(T 3 sg + 1.2 x 10~ 2 s + 1 
W 6 (0, s ) ~ -9.0 x lO^s 3 , + 2.9 x 10~ 2 s + 1 



so 
L5 

so 
1.3 

so 
2.6 

so 
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Clearly these moments converge only when the threshold parameter is very large, so > 10 GeV 2 . Such a large sq is 
irrelevant for X(1600). 
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